Rain initiation time in turbulent warm clouds 
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ABSTRACT 

We present a mean-field model that describes droplet growth due to condensation and collisions 
and droplet loss due to fallout. The model allows for an effective numerical simulation. We study how 
the rain initiation time depends on different parameters. We also present a simple model that allows 
one to estimate the rain initiation time for turbulent clouds with an inhomogeneous concentration of 
cloud condensation nuclei. In particular, we show that over-seeding even a part of a cloud by small 
hygroscopic nuclei one can substantially delay the onset of precipitation. 



1. Introduction 

Droplets start growing by vapor condensation on 
cloud condensation nuclei (CCN) which are typically 
submicron-size particles. In warm clouds, coalescence 
due to collisions contributes the growth until the rain- 
drops (generally exceeding millimeters) fall out of the 
cloud, see e.g. Pruppacher and Klett (1997). Those 
processes can be modelled by the equation for the local 
distribution of droplets over sizes, n(a, t, r) = n{a) and 
the mass content of the water vapor, M(t,r): 



dn(a) . ksM d n(a) 

+ (v • V)n = — - 



+ / da' 



dt po da a 

K(a',a")n(a')n(a") 



(1) 



K(a',a)n(a')n(a) 

2 (or jay 

+ (v • V)M - kAM = -Attspqk f an(a)da . (2) 
dt J 

Here a is the droplet radius, t is time and r is the 
coordinate in space. The first term in the right-hand 
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side of (Jl is due to condensation which (for not very 
large droplets) changes the droplet size a according to 
da 2 1 dt = ksM/pq where k is vapor diffusivity, s is the 
degree of supersaturation and p = 10 3 kg • m~ 3 is the 
density of liquid water. The time needed for condensa- 
tion to grow a droplet of size a is generally proportional 
to a 2 . The second term in the rhs of describes co- 
alescence due to collisions, here a" = (a 3 — a' 3 ) 1 / 3 is 
the size of the droplet that produces the droplet of 
size a upon coalescence with the droplet of size a'. 
The collision kernel is the product of the target area 
and the relative velocity of droplets upon the contact: 
K(a, a')~ir(a+a') 2 Av. According to the recent precise 
measurements (Beard et al, 2002) the coalescence effi- 
ciency of the droplets in the relevant intervals is likely to 
be greater than 0.95 we put it unity in our calculations. 
Collisions change the concentration n(a) on a timescale 
of order 1/ K (a, ai)n(ax) where collision kernel if is a 
fast growing function of droplet sizes. 

Since condensation slows down and coalescence ac- 
celerates as the size of droplets grow then one can in- 
troduce a crossover scale a*, determined by K{a Sf )n ~ 
ksM/ poa 2 . The growth up to the crossover scale is 
mainly due to condensation while coalescence provides 
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for the further growth. The crossover scale a* depends 
on n and is typically in tens of microns (see below). 

To describe the six unknown functions, n,M, s,v, 
one must also add the equation that describes the tem- 
perature change (that determines s) and the Navier- 
Stokes equation for the velocity. Such system cannot be 
possible solved numerically with any meaningful resolu- 
tion neither presently nor in a foreseeable future. The 
main problem is a very complicated spatial structure of 
the fields involved particularly due to cloud turbulence. 
Our aim in this paper is to formulate some mean-field 
model which does not contain spatial arguments at all. 
The requirements to this model is that it must give 
the correct qualitative relations between the parameters 
and reasonable quantitative description (at least within 
the order of magnitude) of the real- world timescales. 
According to the two basic phenomena involved (con- 
densation and collisions), the main problems in space- 
averaging the equations are related to the proper de- 
scription of the two phenomena: mixing and diffusion 
of water vapor and the influence of cloud turbulence 
on collisions. We address them in Sections 0] and [3] re- 
spectively. We use the model to study the evolution of 
n(a, t) starting from sub-micron sizes all the way to the 
moment when droplet fallout significantly decreases the 
water content in the cloud. We shall call this moment 
the rain initiation time and we study how that time de- 
pends on initial vapor content and CCN concentration 
and on the level of air turbulence. 

2. Growth by gravitational collisions 

For the parameters typical for warm precipitating 



clouds {sM/po = 1(T 8 1(T 9 and n = 10 6 



10 a 



collisions are negligible for micron-size droplets (Prup- 
pacher and Klett 1997). For droplets larger than couple 
of microns, Brownian motion can be neglected and the 
collision kernel in a still air is due to gravitational set- 
tling: 



K g {a, a ) =Tr(a-\-a') 2 E(a : a')\u g (a) — u g (a')\ 



(3) 



The fall velocity u g is obtained from the balance of 
gravity force 4irgpoa 3 /3 and the friction F{u gi a). The 
friction force depends on the Reynolds number of the 
flow around the droplet, Re a = u g a/v. When Re a 
is of order unity or less, F — &irvpau g and u g = gr 
where p is the air density and r = {2/§)(p / p){a 2 / v) 
is called Stokes time. We use u g — gr for a < 40^m 
and take u g (a) from the measurements of Gunn and 
Kinzer (1949) for a > 50pm with a smooth interpola- 
tion for 40pm < a < 50pm as shown in Fig^ The 
dotted straight lines have slopes 2, 1, 1/2. One can see 
that u g oc a 2 at a < 40 pm. There is an intermediate 
interval with an approximately linear law u g oc a for 
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Fig. 1. Terminal fall velocity. 



40 pm < a < 400/im. When Re a ;§> 1 one may expect 
F oc pa 2 u 2 as long as droplet remains spherical; that 

gives u g oc y/ agpo/p. Square-root law can be distin- 
guished between 400 pm and 1 mm while the growth of 
u g {a) saturates at larger a due to shape distortions. Hy- 
drodynamic interaction between approaching droplets 
is accounted in K g by the collision efficiency E, which 
values we take from Pinsky et al (2001) at the 750 mbar 
altitude. 

It is of practical use to be able to predict the time 
left before rain starts given the knowledge of droplet 
distribution at a given instant. Such distributions can 
be measured with high accuracy by optical and other 
methods. Drop size distributions measured in many dif- 
ferent types of clouds under a variety of meteorological 
conditions often exhibit a characteristic shape (Prup- 
pacher and Klett, 1997). Generally the concentration 
rises sharply from low value maximum, and then de- 
creases gently toward larger sizes, causing the distri- 
bution to be positively skewed with a long tail toward 
the larger sizes. We approximate such a shape with 
half-Gaussian 6(a — clq) exp(— (a — ao) 2 /2a 2 ) where 6 is 
a step function. We thus characterize the distribution 
by two parameters: the mean position P and the width 
a. Since we mainly consider narrow initial distributions 
(tr < P), the rain initiation time does not depend sub- 
stantially on the initial shape. We start from purely 
gravitational collisions in a still air that is solve the 
space-homogenous version of Q with no condensation 
term: 



dn(a) 



-n(a)- 



(4) 
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K g (a',a")n(a')n(a") 
2(a"/a) 2 



K g (a , a)n(a')n(a) 



The first term in the rhs of models the loss of « 
droplets falling with the settling velocity u g from the 1 
cloud of the vertical size L. Since L are generally | 
very large (from hundreds meters to kilometers) and | 
u g (a) grows with a (see Fig. Q] below) , fallout is rele- | 
vant only for sufficiently large drops (called raindrops) 
with sizes millimeter or more. The collision (Smolu- 
chowsky) term describes the propagation of distribu- 
tion towards large sizes. The asymptotic law of propa- 
gation depends on the scaling of K g (a,a'). If the colli- 
sion kernel is a homogeneous function of degree a [that 
is K g (^a,^a') = £ a K g (a,a')] one can show that for 
a larger/smaller than three the propagation is accel- 
erating/decelerating while for a = 3 it is exponential 
In a oc t (see van Dongen and Ernst 1988; Zakharov et 
al 1991). Our numerics show, however, that the inter- 
vals of sizes a where a is approximately a constant are 
too short for definite self-similarity of the propagation 
to form both for narrow and wide initial distributions. 
This is due to complexity of both functions, u g (a) and 
E(a,a'). We thus focus on the most salient feature of 
the propagation, namely study how the amount of wa- 
ter left in the cloud, W, depends on time. The decrease 
of that amount is due to a concerted action of collisions 
producing large drops and fallout. 

The droplets radii space was discretized, i.e. the 
droplets size distribution n(a, t) was presented as the set 
of concentrations rn(t) of droplets with radius oj. The 
total mass of vapor and water in droplets is conserved 
in our calculations. The grid of radii was taken approx- 
imately exponential at sizes that are much larger than 
the size of initial condensation nuclei, with 256 points in 
unit interval of natural logarithm. The collision term in 
Smoluchowsky equation was treated as follows: let the 
radius (af -l-a:-) 1 / 3 of the droplet resulted from merging 
of the two with radii di and a,j to be in between of two 
radii au and a k +x from the grid. Then the collision re- 
sults in decreasing of rii and nj by quantity dN that is 
determined by the collision kernel, while the concentra- 
tions iik and rik+i are increased in such a way that sum 
of their change is dN and the whole amount of water in 
droplets is conserved in coalescence: 

Srii — Srij = —dN = —5rik — 6nk+i , 



a\6n k + a\ +1 5n k+1 = (a 3 + a 3 )aW . 
Sn k+ i — dN(af + a 



4)/(a 3 

6n k = dN(4 +1 - a 3 + a))/{a k 



fc+i 

3\ /f„3 



(5) 



The total amount of water (the sum of the part that 
left and remained in the cloud) was a conserved quan- 
tity, up to 10 -6 accuracy, during the whole simulation. 




Fig. 2. Fraction of water left in the cloud as a function 
of time. P — 13 /im. 



Note that our scheme automatically keeps the numbers 
positive: if dN is greater than either m or rij, then we 
choose dN = min{ni, rij}, so that rij and rij are also 
not negative after every elementary collision process. 
Let us stress that our scheme is conservative both in 
mass and number of droplets (comparing to the non- 
conservative scheme of Berry and Reinhardt, 1974 and 
the scheme of Bott, 1998 which was conservative only in 
mass and needed special choice of the time step to keep 
positivity). The minimal time step needed for our cal- 
culations was estimated from characteristic timescales 
of our problem to be 0.1 s. We have checked that the 
decrease of the time step below dt = 0.05 s does not 
change the results, the figures below all correspond to 
that dt. 

The graphs W(t) are shown at FigsEJand|SI(for L = 2 

km) and they are qualitatively the same both for narrow 
and wide initial distributions. At the initial stage, W 
decreases slowly due to the loss of drizzle. After large 
raindrops appear, loss accelerates. At every curve, the 
star marks the moment when respective d 2 W/dt 2 are 
maximal. After that moment, the cloud looses water 
fast so it is natural to take £* as the beginning of rain. 
Figure 0] shows how the mass distribution over sizes, 
m(a) oc a 3 n(a) evolves with time. One can see the ap- 
pearance of secondary peaks and distribution propagat- 
ing to large a. The moment t* seems to correspond to 
the highest value of the envelope of the curves m(a, t) 
of the coalescence-produced drops. One can see from 
Figure 01 that the peak at mass distribution is around 
200 microns and most of the droplets are below 500 mi- 
crons at t = i*. The same character of the evolution 
W(t) can be seen in the next section for the ab initio 
simulations of 11121) . 

The rain initiation time t* defined in that way is 
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Fig. 3. Fraction of water left in the cloud as a function 
of time. P = 16 /im. 




Fig. 4. Mass density of water, t* ~ 1500 s. 
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Fig. 5. Rain initiation time as function of the width of 
initial distribution a for different initial positions P. 

presented in Figures [5] and [5] against the width and the 
position of the initial distribution. Note the dramatic 
increase in t» with decreasing a for P = 13 fim. The 
mean droplet size P = 14 /im is sometimes empirically 
introduced as the minimal size required for the onset of 
precipitation (Rosenfeld and Gutman 1994). Figures 
and support that observation, they indeed show that 
i* grows fast when P decreases below that size but only 
for very narrow initial distributions and of course there 
is no clear-cut threshold as t*(P) is a smooth (though 
steep) function. The timescales (from tens of minutes to 
hours) are in agreement with the data obtained before 
(see Pruppacher and Klett, 1997, Chapter 15; and Se- 
infeld, J. and S.Pandis, 1998, Chapter 15 and the refer- 
ences therein). FigureElalso shows that for 15 /im P, 
the function t*(P) can be well-approximated by a power 
law oc P^ 1 with 7 w 3. The rain initiation time de- 
pends on the cloud vertical size almost logarithmically 
as shown in Fig. \7\ we do not have an explanation for 
this functional form. 

Here we treated the position and the width of the 
distribution as given at the beginning of the collision 
stage. But of course the distribution is itself a product 
of condensation stage so we now turn to the considera- 
tion of the full condensation-collision model. 

3. Condensation and collisions 

We consider now the space-homogeneous system 



dn(a) ksM d n(a) , ,U a (a) 

n[a) ^ 



dt- 



ps da a 



(6) 
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Fig. 6. Rain initiation time as function of the position of 
the initial distribution. 
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Fig. 7. Rain initiation time as function of the cloud ver- 
tical size. 



If one substitutes here gravitational and Brownian col- 
lision kernels (taken, e.g. from Pruppacher and Klett 
1997) and start from n = 10 7 -j- 10 s m -3 sub-micron 
droplets growing in a medium with sM / p$ = 10~ 8 -j- 
10~ 9 then H6I7|) give unrealistically large rain initiation 
time. The reason for that is well-known: during the 
condensation stage, the distribution shifts to larger sizes 
while keeping its small width over a 2 . For narrow dis- 
tributions, gravitational collisions are suppressed (since 
all droplets fall with close velocities) as we have seen in 
the previous section. Collisions of droplets with similar 
sizes are provided by an inhomogeneous air flow. The 
velocity gradient A provides for the kernel K s — X(a+a') 3 
derived in Saffman and Turner (1956). However, typical 
velocity gradients in the air turbulence (A ~ 10 — 30 s _1 ) 
also do not provide enough collisions (see e.g. Prup- 
pacher and Klett 1997; Jonas 1996; Vaillancourt and 
Yau 2000, and the referenced therein) . Regular ver- 
tical inhomogeneity of supersaturation due to temper- 
ature profile does not broaden n(a) much even with 
the account of turbulence-induced random fluctuations 
(Korolev 1995; Turitsyn 2003). Spatial inhomogeneities 
in vapor content M due to mixing of humid and dry 
air still remains a controversial subject (see. e.g Prup- 
pacher and Klett 1997; Baker et al, 1980) and proba- 
bly can be neglected in cloud cores. We address the 
turbulent mixing of vapor in Section ^considering par- 
tially seeded clouds. We address the turbulent mix- 
ing of vapour in Section 21 considering partially seeded 
clouds. As far as collisions are concerned, the main ef- 
fect of spatial inhomogeneities seems to be the effect of 
preferential concentration that is of turbulence-induced 
fluctuations in droplets concentration (see Maxey 1987; 
Squires and Eaton 1991; Sundaram and Collins 1997; 
Reade and Collins 2000; Grits et al 2000; Shaw et 
al 1998; Kostinski and Shaw 2001, and the references 
therein). We use here the results of the recent the- 
ory by Falkovich et al (2002). Namely, we multiply the 
Saffman- Turner collision kernel K s by the enhancement 
factor (n 1 n2) / (ni) (n 2 ) (which accounts for the effects 
of inertia and gravity in a turbulent flow) and add the 
collision kernel due to the so-called sling effect (droplets 
shot out of air vortices with too high centrifugal accel- 
eration), see Falkovich et al (2002) and Falkovich and 
Pumir (2004) for the details. The total collision kernel 
due to turbulence normalized by the homogeneous ex- 
pression 8Aa 3 factor is presented in Fig.[S]for Re = 10 6 . 
We see that role of concentration inhomogeneities can 
be substantial in the interval between 30 and 60 /j,m. 
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Fig. 8. Turbulence collision kernel normalized by 8Aa 3 
for equal-size droplets at Re = 10 6 . From bottom to top, 
A = 10, 15 and 20 s" 1 . 



The role of sling effect is not significant at those levels 
of turbulence: for the upper curve it gives the contri- 
bution of order of 10% between 25 and 35 /im. 

The system I|6I7I) is our mean-field model where the 
only memory of spatial inhomogeneities are the fallout 
term and the renormalization of the collision kernel K. 
As we show here, this model gives the rain initiation 
times with reasonable quantitative values and proper 
qualitative behavior upon the change of parameters. 
Let us discuss first how depends on n. Here, the 
most important feature is the existence of the mini- 
mum in the function i»(n). That can be explained by 
the competition between condensation and collisions. 
Increasing n one decreases a* (the crossover size for 
which condensation time is comparable to the time of 
collisional growth) and thus decreases the time needed 
for droplet growth. This works until a» is getting com- 
parable to a c ~ (M/npo) 1 ^ 3 - Indeed, when droplets 
grow comparable to a c vapor depletion slows and then 
stops condensation. If one takes the initial concentra- 
tion even larger so that a c < a* then vapor depletion 
stops condensation earlier and collisions are slower for 
droplets of the smaller size a c so that the overall time 
of droplet growth is getting larger. The concentration 
n* that corresponds to the minimal time can be found 
from the (implicit) relation a c ~ a* which corresponds 
to 



(M/n^o)- 1 /^ (M/n^ ) 1/3 ) 



(8) 



That tells that 



M and if K 



a a then 

n* oc s 3 /( 1_a ). One can argue that for small concen- 
trations (generally for maritime clouds), n < n», times 
of condensation and collision stages are comparable. 
Therefore, is a function of the product Ms. For a 
homogeneous kernel, i* oc n~ 2 ^ 2+a ^> (Al s)~ a ^ 2+a \ For 



large concentrations (generally for continental clouds), 
n > n*, the rain initiation time is mainly determined 
by collisions so it is getting independent of the super- 
saturation and oc n ^ a ~' 3 ^ 3 M~' a / 3 . 

By numerically solving Hfil7[) with both gravity and 
turbulence-induced collisions we obtain the rain initia- 
tion time (also defined by the maximum of cPW/dt 2 ) 
as a function of the CCN concentration n for different 
values of the supersaturation s and the vapor content 
M. The grid of radii was approximately exponential at 
sizes that are much larger than the size of initial con- 
densation nuclei (with 200 points in unit interval of nat- 
ural logarithm) . The condensation of vapor was taken 
into account by working on evolving grid of radii a,i(t) 
keeping conserved the total mass of water in droplets 
and vapor. Collisions were modelled according to (0) 
described above. Note that the numerical scheme we 
employ here has an additional advantage (comparing to 
those described in Pruppacher and Klett, 1997; Berry 
and Rcinhardt, 1974; Bott, 1998) of accounting simulta- 
neously for condensation and collisions while respecting 
conservation laws. We used the time step dt = 0.01 s 
during the condensation phase, on a later stage (domi- 
nated by coalescence) dt = 0.1 s was enough. Those re- 
sults are presented in FigureElfor L = 2 km and A = 20 
s . The solitary point at the lower part corresponds to 
M = 6 g- m~ 3 , s = 1/150. The three solid lines corre- 
spond to AI = 3 g- m~ 3 while the three dashed lines to 
M = 1.5 g- m~ 3 . Inside the triplets, the lines differ by 
the values of the supersaturation, from bottom to top, 
s = 1/75, 1/150, 1/300. We see that indeed the graphs 
i* (n) all have minima. The position of the minimum is 
proportional to M as expected and approximately pro- 
portional to s" 1 / 2 which would correspond to a ~ 7 
in this interval of sizes. We see that the left parts of 
different curves with the same product sAI approach 
each other as n decreases. To the right of the minima, 
the curves with different s but the same M approach 
each other as n increases. That supports the previous 
conclusions on the respective roles of condensation and 
collisions in determining the rain initiation time. 



4. Delaying rain by hygroscopic over-seeding 

That the rain time is a non-monotonic function of the 
concentration of droplets may provide a partial explana- 
tion for the conflicting observations of the effect of hy- 
groscopic seeding. By seeding clouds with hygroscopic 
aerosol particles one can vary the number of cloud con- 
densation nuclei and thus the number of small droplets 
at the beginning of the cloud formation. It was ob- 
served that such seeding in some cases suppresses pre- 
cipitation (see e.g. Rosenfeld et al 2001), while in other 
cases enhances and accelerates it (Cotton and Pielke, 
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Fig. 9. Rain initiation time as function of CCN concen- 
tration n for different supersaturations s and vapor contents 
M. 

1995; Mather 1991), see also Bruintjes (1999) for a re- 
cent review. 

It is often desirable to postpone rain, for instance, to 
bring precipitation inland from the sea. The fact that 
i* grows when n > n* suggests the idea of over-seeding 
to delay rain. This is considered to be unpractical: "It 
would be necessary to treat all portions of a target cloud 
because, once precipitation appeared anywhere in it, 
the raindrops . . . would be circulated throughout the 
cloud ... by turbulence" (Dennis, 1980). We think that 
this conclusion ignores another, positive, aspect of cloud 
turbulence namely the mixing and homogenization of 
partially seeded cloud during the condensation stage. 
Let us describe briefly how it works for two cases. 

Consider first seeding a part of the cloud comparable 
to its size L c . Note that we do not consider here adding 
ultra-giant nuclei, we assume seeded CCN to be com- 
parable in size to those naturally present. According 
to the Richardson law, the squared distance between 
two fluid parcels grows as et 3 so that the rms difference 
of vapor concentrations between seeded and unseeded 
parts decreases as t~ 9 ^ 4 when t 3 > i 3 , — L 2 /e (e is 
the energy dissipation rate in turbulence). To see how 
different rates of condensation interplay with turbulent 
mixing we generalize the mean-field system Il(il7t de- 
scribing seeded and unseeded parts by their respective 
7ii,ri2 and x\ — s\Mi, X2 — S2M2 and link them by 
adding the term that models the decay of the differ- 
ence: dxi/dt = ... — (xi — xj)t(t + t )~ 2 (9/4). As a 
crude model, we assume two parts to evolve separately 
until t = 2to, then we treat the cloud as well- mixed and 
allow for the collisions between droplets from different 
parts. That actually underestimates the effect of seed- 
ing and can be considered as giving the lower bound 
for the time before rain. The results of simulations are 
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Fig. 10. Fraction of water left in the cloud as a function of 
time (A). Time of rain t, as a function of CCN concentration 
no (B). The lower part (boxes) corresponds to a half-seeded 
cloud (the half-sum of concentrations is used as abscissa). 



shown in Fig. [H] for to = 180 s, L = 2 km and A = 20 
s . It is seen from Fig. I10t \ that the water content 
W changes similarly to what was shown in Figs. I2l.il 
and the rain initiation time is again determined by the 
maximum of d 2 W/dt 2 . The respective times are shown 
against n = (n 1 +n 2 )/2 by boxes in Fig.llOb. The time 
increase is less than for homogeneous seeding but is still 
substantial. The fraction of the cloud still unmixed af- 
ter the time t decreases by the Poisson law exp(— t/to). 
Taking n\ = 100 cm -3 one sees that for a time delay of 
10 min one needs to seed by n-i ~ 3000 cm~ 3 . 

Second, consider seeding by N particles a small part 
of the cloud which (unseeded) had some no and would 
rain after £*. After time i* the seeds spread into the 
area of size {et 3 ) 1 / 2 with the concentration inside the 
mixed region decaying as n(i*) = N{et 3 )~ 3 / 2 (for strat- 
iform clouds one gets N(et 3 )^ 1 ). To have an effect of 
seeding, one needs n(t*) > no which requires N > 10 15 
for no = 50 cm -3 , f* = 10 min and e = 10 cm s , 
With sub-micron particles weighing 10~ n g that would 
mean hundreds of kilograms which is still practical. 
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5. Summary 

We believe that our main result is a simple mean-field 
model H6I7|I which demonstrates non-monotonic depen- 
dence of the rain initiation time on CCN concentration. 
As the CCN concentration increases, the rain initia- 
tion time first decreases and then grows as shown in 
Figs. I9lf 01 The simple modification of this model for 
an inhomogeneous case described in Sect. 0] shows that 
one can increase the rain initiation time even for a cloud 
partially seeded by hygroscopic aerosols. 

We acknowledge support by the Ellentuck fund, by 
the Minerva and Israel Science Foundations and by NSF 
under agreement No. DMS-9729992. We are grateful 
to A. Khain, M. Pinsky and D. Rosenfeld for useful 
discussions. 
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